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We analyze the frequency-dependent noise of a current through a quantum dot which is coupled to 
Fermi leads and which is in the Coulomb blockade regime. We show that the asymmetric shot noise 
as function of frequency shows steps and becomes super-Poissonian. This provides experimental 
access to the quantum fluctuations of the current. We present an exact calculation for a single dot 
level and a perturbative evaluation of the noise in Born approximation (sequential tunneling regime 
but without Markov approximation) for the general case of many levels with charging interaction. 



The shot noise is a striking consequence of charge 
quantization and allows to characterize the transport 
of individual electrons The symmetry of the noise 
S{lu) is important: For a classical stationary system, the 
noise (for autocorrelations) is always symmetric in the 
frequency lo. However, for a quantum system, the noise 
can by asymmetric due to the non-commutativity of cur- 
rent operators at different times. It was recently found 
that such an asymmetric noise can be detected since the 
noise frequency lo corresponds to a quantum of energy fku 
which is transferred from the system to the measurement 
apparatus |^ H, Q which was demonstrated experimen- 
tally n. This means that antisymmetric quantum effects 
in the noise can be measured and isolated from the clas- 
sical (symmetric) effects. We show that striking asym- 
metric effects appear in the shot noise S{uj) of a quantum 
dot with steps as function of co, giving a supcr-Poissonian 
Fano factor. Our analysis is based on a perturbative ap- 
proach which remains valid in the quantum limit with 
large lo (where a Markov approximation typically invoked 
would not be valid). We confirm our perturbative results 
by an exact calculation of the noise for a dot with a single 
level. We note that quantum dots are good candidates 
for an experimental test of our predictions since such sys- 
tems have been studied extensively over the years, both 
experimentally and theoretically [a la S Hy? 113 • 

We consider the operator /; which describes the current 
in a lead We define the current noise. 
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in terms of the (non-symmetrized) correlation function, 
{Ii{t)Ii') — Ti Ii{t)Iii p. Here, p is the stationary density 
matrix (of the full quantum system). Note that Ii{t)Ii' is 
not a Hermitian operator and thus does not correspond to 
a classical observable. How should we interpret Eq. ^7 
On the one hand, one can avoid the non-Hermitian op- 
erator by arguing heuristically that Eq. (QJ is "unphysi- 
cal" and by instead considering the correlation function 
in terms of the symmetrized operator i [/;(t)/;/ + /;///(t)] 
p^ . One then obtains the symmetrized noise, S'f;^™(w) = 
^[Sii'{uj) + Si'i{—uj)]. For I = V , this noise does not 
depend on the sign of w. This corresponds to a mea- 
surement apparatus which does not discriminate 
between absorption and emission of energy by the sys- 



tem, and thus cannot detect all noise characteristics. On 
the other hand, we have (//(—<)/;') = {Ii'{t)Ii)* and thus 
Sw{lu) = Sinluj)*, since /; is Hermitian and p stationary. 
Thus, Su (lo) is a real quantity which can be regarded as 
an observable. Indeed, there are setups, where Su{uj) can 
be accessed experimentally |2,[3[0,|3- This justifies to 
consider the asymmetric quantum noise [Eq. (^] as we 
do in the following. 

Quantum dots. To illustrate the presence of asymmet- 
ric shot noise contributions due to quantum effects, we 
consider now a concrete system of a quantum dot in the 
Coulomb blockade regime coupled to Fermi leads 
I = 1,2,... at chemical potentials pi. When only a sin- 
gle dot level is present, the noise can be calculated even 
exactly (see below). This is however not possible for 
systems with many levels and charging interaction, for 
which we now develop a perturbative approach. We as- 
sume weak coupling such that current and noise are dom- 
inated by the sequential tunneling (ST) contributions, 
valid for kT > 7 with temperature T and level width 
7. We model the combined system with the Hamiltonian 
H = i^ioad + + Ht, which describes leads, dot, and 
the tunnel coupling between leads and dot, resp., and 

with Ho = i/icad + Hd. We let i?ioad = J2lkcr <^lkc\f,„Clka , 

where c\j^^ creates an electron in lead I with orbital state 
k, spin cr, and energy ejfc. The electronic dot states 
\n) are described by -ffdl'^-) = En\n), including charging 
and interaction energies. We use the standard tunneling 
Hamiltonian Ht = I];pfc<j ^Ip^.ka'^pc' + with tunnel- 
ing amplitude and where d'^piy creates an electron on 
the dot with orbital state p and spin a. The state of 
the combined system is given by the full density matrix 
p, while the electronic states of the dot are described by 
the reduced density matrix, ps = TrR p, where the trace 
is taken over the leads. We assume that at some initial 
time to the full density matrix factorizcs, p(io) — p'sPb.^ 
with the equilibrium density matrix of the leads, pp. 
From the von Neumann equation p = —i [H, p] one finds 
[1^ the generalized master equation for the reduced den- 
sity matrix, psit) = -tLsPs{t) - dt' M{t')ps{t-t'), 
where the kernel M is the self-energy superoperator. 
Since we consider the weak coupling regime, we pro- 
ceed with a systematic lowest-order expansion in Ht- 
We obtain M(t) = Tr^ Lt e~^^°^ LtP^, where we define 
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the superoperators LsX = [Hd, X], LqX = [Hq, X], 
and LxX = [Ht, X]. In the foUowing, we work m 
the Laplace space, f{t) t-^ f{u}) = dt e*"*/(t) (we 
take Imw > but our results remain well-defined for 
Imw — > 0). Then, the time evolution of ps reads 



ps{h) - iujps{(^) = M{uj) ps{(^) 



(2) 



with A4{u!) = —iLg — M{uj) and with the lower boundary 
of the Laplace transformation shifted to Iq. We take to 
—oo and assume that the system has relaxed at the much 
later time t = into its stationary state pg = ps{Q) = 
\\m^^^Q{—iw)ps{oj). We multiply Eq. (01 by — icj, take 
^ 0, and find the equation Ai{Q)ps = from which 
we get Ps- 

Current. We calculate the current /( flowing from 
the dot into lead I and vice versa. The current oper- 
ators are Ii{t) = (— l)'eg;(t) = (— l)'ie[iJT, <?;(*)] where 
qi ~ ^^.g. c\f.^cika is the number of electrons in lead We 
choose the sign of /; such that (/i) = {I2) in the case of 
two leads. We now introduce the projectors P = p^TrR 
and g = 1 - P with the properties PLtP = Q ^ PIiP 
and PLq = L^P. We evaluate (/;) by inserting P -\- 
g = 1 and find {hit)) = Tr /zge-^^^o+^^^^'-'oVl^o) = 
-iTr/i/j*^ dt' Qe-'^°'^*-*">LTPp{t') + 0{H^). This mo- 
tivates introducing the following superoperators, = 



W> 



with W> {t) 



and W<(t) = -iTTnLre-'^^^Iipl, and I^m'(t) 
TrR 6"'^""^/;/^^ . Note that these superoperators act 
only on the dot space, which considerably simplifies fur- 
ther evaluations. In the ST regime, the current is 



Here, uj is real and the limit w ^ is well behaved [the 
S{uj) contribution is cancelled by (/;)(/;') in Eq. For 
a superoperator S, we have defined S such that {SA)^ = 
SA'', thus Sbcinm ■= {b\{S\n){m\)\c) = (5cb|mn)*, and 
S{lu)A = J^dte"^'^[S{t)A]K For deriving Eqs. ®-l|3J), 
we have made no Markov approximation where we would 
evaluate M(i)e'*^° at w = and equivalently for the 
other superoperators We note that for the current 
[Eq. @], M and are both evaluated at w = 0, thus 
the stationary current does not contain non-Markovian 
effects. However, the Markov approximation changes the 
noise, which for I = I' becomes symmetrized, Sfi^^°^{uj) ~ 
^Mkov^_i^-j particular, the antisymmetric noise 

contribution (showing pure quantum effects) cannot be 
obtained in the Markov approximation. 

We now return to the exact expression of the noise in 
Born approximation [Eq. I^J] and explicitly calculate the 
matrix elements of the various superoperators, 

-M(t)e'*^V5 = E i^+Ps - 9+ Ps) + h.c, (5) 



W<>\t)e^*^« = (-l)'e( T £), 



W^,{t) 



Sw 9+, 



(6) 
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with G?t = G''°"*(i) ± G'''"(t) and gi^ = 
j:bn\b){n\TTs{G'±\n){b\} M We define 

= V27rt';<j Z]p^i^("-MpT|TO), with spin depen- 
dent density of states vi^ in lead /. The matrix elements 
of the remaining superoperators are 



(Ii) ^TvsW!{uj = 0)ps. 



(3) 



This implies that the superoperator VF/ corresponds to 
the current. 

Quantum shot noise. We evaluate the noise [Eq. 
^] in lowest order in Ht but without any further 
approximation. It is sufficient to consider t > 0, 
since {Ii{-t)Ii') = {Ii>{t)Ii)*. Using again P + 
g = fl we get {Ii{t)Ii,) = Tr/,ge-*^*P/rgp + 
Tt IiQe~''^*QIi'P p, while we neglect the higher-order 
term IiQe'^^^QIvQ p in Ht- The goal is to fac- 
tor out one of the following expressions. First, we ex 



M 0;(A6„-c.) + ^), (8) 



with A6„ ^Eb-En, and p± = log{27rfcT/[(l T l)ec/2 ± 
Abn-uj]} + Rei;[^+i (Ab„ Tto-fn) /2nkT] . Here, ijj 
is the digamma function. The terms arise from 
the principal values V defi{e)/{€ — Ah„ + lu) and 
r J^" de [1 - /;(e)] /(A„6 -e + u) with bandwidth cutoff 
pand Qe-'^Q'Q = Qe'^^^^Q + 0{Ht) to leading or- If neglect uj with respect to the large energies A^^ 

der. Second, we evaluate the conditional time evolution and - A„b, the first term of p (and thus e^) drops out 
p^sit) := TrRe"*-^VR- We find that p| is the formal 
solution of Eq. O) with initial value 1, thus p'si^) = 
-[iw + 7W(w)]-i Jli. Finally, we use TrR /,g g-'-^Vr = 
J* dt' W> (t - t') p's it') and TrR/,ge-'^*g/i,pO = 

wi,i' it) + /„* dt' j^'dt" w> (t tn it' n w< it"). 

We obtain the noise correlation in the ST regime, 

Su'iuj) = 2Trs [Wliu) p%iu) [VF>(0) -t- W< iuo)] + 
+ Wli-u;) p%i~u) [#>(0) + W< i-uj) 
+ Wi,viuj) + Wv,ii-uo)^ps. (4) 



in Sii'iuj). We note that the contribution corresponding 
to Eq. ^ has been calculated for the symmetrized noise 
of a single electron transistor with a continuous spectrum, 
using a phenomenological Langevin approach [13 . With 
these results, Eqs. and jSJl-lO, it is straightforward 
to find Sw iuj) for an arbitrary dot spectrum; one only 
needs to evaluate simple algebraic expressions. 

We now identify the regime where the asymmetric 
noise properties become most apparent. Asymmetries 
arise from the w dependence of Eqs. ((SI)-®, i.e., are 
most prominent for \uj\ > kT, with steps occurring at 
\uj\ ~ |Af„j — fill (see below). In this regime, the Markov 
approximation breaks down [l^ and the noise probes 
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Figure 1: Quantum dot coupled to two leads and in the 
sequential tunneling regime, (a) Large Zeeman splitting, 
Az > A/i + u. When the dot is empty and lu > — fi2, 
an electron from lead 2 absorbs energy uj and tunnels for a 
short time onto the dot, contributing to the noise 822- Sim- 
ilarly, for u > fj,i — El, the electron on the dot can tunnel 
into lead 1, contributing to 5*11. (b) Smaller Zeeman split- 
ting, here Az = A/i/4. When the energies Ef^i ± u} (dotted 
lines) are aligned with /ii,2, the shot noise has a step, see Fig. 
|l(solid line). 



non-Mar kovian effects. Further, since the ST regime 
is valid for kT > \t\'^, we have w > Thus, the 

conditional time evolution becomes P%{lli) ~ i/u, since 
M{lu) ~ t'^. For the noise [Eq. only the two last 
terms are relevant, since they are of order while the 
other terms are of order t'^ /ui and can be neglected. Less 
formally, this is because multiple tunneling processes [de- 
scribed by pg{w)] do not occur on the short time scales 
corresponding to large tu. Thus, only the individual (un- 
correlated) tunneling events contribute, leading to shot 
noise. 

Next, we discuss specific cases, see Fig. where a dot 
is coupled to two leads / = 1, 2 and a voltage bias = 
fii — /i2 is applied. We assume single energy level spacing 
and Coulomb charging energy larger than temperature, 
bias, and noise frequency. We consider the dot states |0), 
with an even number of electrons and state jcr) where 
an electron with spin a =t,i is added to the dot [20l| . 
For an applied magnetic field B, the Zeeman splitting is 
Az = gpBB = El— E-^ > 0. We consider the ST regime, 
pi > E^ > /i2, and define the tunnehng rates 7^ = l^g'^p. 

Dot with single level. First, we assume a large Zee- 
man splitting such that only the spin ground state 1 1) 
is relevant, see Fig. Q^a), and we omit the index "f. 
Since in this regime only one dot level is involved, there 
are no charging effects between different levels. Thus, 
i?d = E^d'^d and so the full Hamiltonian H is bilin- 
ear and can be solved exactly. The symmetrized noise 
was calculated for this system and discussed for lo — Q 
0. We now calculate the asymmetric noise [Eq. J^l] 
for finite w exactly. For this, we solve the Heisen- 
berg equations for d{t) and Qfc(i) and find the cur- 
rent operator, /i(t)/e(-l)' = Y.kVk' h'^lk'^ik + H.c. + 

Y.i'k'i"k"ilil\ti?) jk' i*knc\,yCink" . Here, the lead op- 
erators cik are evaluated at time (i.e., when the 
leads are at equilibrium) and we have defined = 
it*ti, e'^'''''''-"'''>'-*-*°'>/ {ei,k' - El -i-f) and 7 = (71 + 
72)72. Now we insert Ii{t) into Eq. ^ and readily obtain 
the asymmetric noise, containing all quantum effects. We 



consider the coherent non-pcrturbativc regime of strong 
coupling to the leads in the quantum limit of large fre- 
quencies, w > 7 > kT. We obtain the shot noise 

^n^'^) = E ^ ± T mO Hpi') ~ Hpi t lo)] , 



(10) 

where h{e) — arctan [(e — £'|)/7]. Note that the noise 
shows steps at uj = ±\Ei — Hi] with width 7 [2l]|. Fur- 
thermore, for uj > \Ei — fj,i\, Afj,, the noise is asymmetric 
and saturates at S^{uj) = e^7;, while S^{~uj) — 0. 

Let us now consider the ST regime kT > 7 in the exact 
solution. For a; > 7, we then find 

Sni^) = E ^ [^i.Ti ± fi'im [Si,±i T fi{Ei ± Lo)] . 
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(11) 

Again, the noise shows a pronounced asymmetry. We 
can now compare Eq. (ini) with the noise obtained in the 
perturbative approximation [Eq. Ql] and find that they 
agree. We further consider cu > Ap + kT in Eq. such 
that MEi+u;) = and fi{Ei - w) = 1, leaving MEi) 
unrestricted. In this case, the (asymmetric) shot noise is 



(12) 



whereas Sfi{—uj) = 5*^2 (±w) = 0; this is the same result 
as we have found for strong coupling [Eq. Hl()(l ]. The 
interpretation is that for Sfi{—Lo) the detector absorbs 
energy which, however, cannot be provided by any 
tunneling process. On the other hand, for Sii{Ld) the 
detector provides energy uj. Thus, if the dot is empty, an 
electron with energy E^ — ld can tunnel from the Fermi 
sea I into the dot, and, if the dot is filled, an electron 
can tunnel from the dot into an unoccupied lead state 
of energy E^ -\- to, see Fig. QJa) In both cases, the 

tunneling occurs with rate 7; and thus the contribution to 
the autocorrelation is 6^7; 6{t). Note that for \Ei — > 
kT, the noise is 5'fi(tj) = e(/) (71 -I- 72)/72- Thus, for 
large uj, the frequency dependent Fano factor, _F'ii(ti;) = 
Sii{uj)/e{I), is 2 for 71 = 72, and can even become larger 
for 71 > 72, in contrast to the Markovian case where we 
find it to be 1. Thus, we find that the quantum shot noise 
is super-Poissonian. Moreover, away from the ST regime, 
say for E^ + kT > fj.i, the dot remains always in state 
|0) and only a small (higher-order in Ht) cotunncling 
current (J) fiows through the dot However, the noise 
can still be of lower order, it is Sfi{uj) = e'^^ifi{Ei — uS) 
for large \lo\, resulting in a large Fano factor Fu(lo) and 
super-Poissonian shot noise. 

Dot with two or more levels. Second, we consider the 
regime where the state 1 1) becomes relevant and charging 
interaction enters (here no exact solution is available). 
We consider a small Zeeman splitting such that fii > 
Ei^l > fi2 and fi{Ei) « fi{Ei), see Fig.UIb). Using Eq. 
Q, we calculate the noise Si^{uj) and plot it in Fig. [21 
(solid line). For large |a;|, such that /;(i?cr + |a;|) = and 
fi{Eo- — \uj\) = 1, the noise vanishes for w < while for 
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Figure 2: The Fano factor Fii{uj) = Sii{uj)/e{I) in the 
shot noise regime Ajj, > kT as function of noise frequency 
LV. (Asymmetric noise at frequencies up to 90 GHz has been 
measured 0.) We consider T = 100 mK, Afi/e = 460 fiV, 
= im + H2)/2, 71 = 72 = 5 X 10^ s-\ and 5 = 2. We 
use the full expression for the noise S\i [Eq. (solid line) 
and within Markov approximation, 5'ii^''°^, (dotted line), for 
B = IT [see Fig.mb)], thus = A/i/4 and (/) = 530 pA. 
We also show Sa (dashed line), being strongly asymmetric, 
where B = 3T, = 3A/i/4, and (/> = 400 pA. The dip 
near a; = is due to the charging effect of the dot, while the 
steps at uji (see text) arise from additional transitions for in- 
creasing uj and provide a striking effect in the quantum shot 
noise. Note that these steps disappear when the noise S\-i is 
symmetrized (some features remain for 71 ^ ^2). 



cij > it saturates at 



71 + 72 



7i[1 + /i(-Bt)]+72[1 + /2(-Bt)]- 



(13) 



More generally, for the weaker assumption 
> 7, the numerator in Eq. H13() becomes 
\ Ei',±,<. 7r [^i,Ti ± h'{E{)\ [<Si,±i T h{Ea ± c^)]. Thus, 
the noise S\-^ shows four steps at uJi = ±(/ii — i?t,x), 
corresponding to the dotted lines in Fig. ^b). The 
interpretation is that for increasing w, more energy is 
available and more tunneling processes are allowed. 
Namely, if w > — (/ii — i5f), an electron with spin "f 
from lead 1 can emit energy lo and tunnel onto the 
dot, and if a; > fii ~ E^-, an electron on the dot with 



spin 1 can absorb energy u and tunnel into lead 1. 
Analogous processes occur for spin |. The steps in 
Sfi are broadened due to and temperature, and the 
step is oc tanh[{w-aj^)/2kT]. For Apt > kT, A:,, 
i.e., fi{Ea) = (5ii, the height of the steps in the Fano 
factor Fii(tj) is \ for uji < Q and 71/272 for oji > 0. 
(These heights change for spin-dependent tunneling, 
7^^ ^ 7j^.) Next, we consider an intermediate Zeeman 
splitting. El > idi + kT. In this regime (c), the dot is 
either in state |0) or while the state 1 1) is never 
occupied and so no additional tunneling process occurs 
for uj > —{El — fii). Thus, the steps in the noise are at 
u! — ±(/ii — El) and a,t uj = Ei — fii, see Fig. |5] (dashed 
line). For large uj, the noise saturates at 



/' = 1,2 



Here, we have allowed for spin-dependent tunneling. If 
we exclude the contributions involving state | ],) by set- 
ting 7;''' = 0, we recover Eq. (|12|) . Generally, we see 
that the shot noise 811(10) of a quantum dot consists of 
a scries of steps and is monotonically increasing, apart 
from features near a> = 0. Each dot level with energy 
(chemical potential of the dot) Ej gives rise to steps 
at ±(/i; — Ej) if the level is inside the bias window, 
/zi -I- kT > Ej > /i2 — kT, and to a single step at — Ej \ 
otherwise. We stress that such a highly asymmetric 
Sii{lu) can be observed with an appropriate measurement 
apparatus For sufficiently large \uj\, the anti- 

symmetric contribution becomes 5 [Su{llj) — Su{—uj)] = 
isign(a;)S';i(|a;|) and is, e.g., given by Eqs. lfTT|l - l(THl . 

In conclusion, we have derived the asymmetric shot 
noise of a quantum dot exactly for a single dot level and 
in the weak coupling regime for many levels. We have 
shown that the shot noise exhibits strong asymmetric and 
super-Poissonian effects in the quantum limit. We thank 
F. Marquardt, W. Belzig, G. Burkard, A. Cottet, J.C. 
Egues, H. Gassman, D. Saraga, and C. Schonenberger 
for discussions. We acknowledge support from the Swiss 
NSF, NCCR Nanoscience Basel, DARPA, and ARO. 
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